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The Bloch sphere is a generic picture describing 
a coupled two-level system and the coherent dy- 
namics of its superposition states under control of 
electromagnetic fields^. It is commonly employed 
to visualise a broad variety of phenomena rang- 
ing from spin ensembles 2 and atoms^ to quantum 
dots 4 and superconducting circuits 5 . The under- 
lying Bloch equations^ describe the state evolu- 
tion of the two-level system and allow character- 
ising both energy and phase relaxation processes 
in a simple yet powerful manner™. 

Here we demonstrate the realisation of a 
nanomechanical two-level system which is driven 
by radio frequency signals. It allows to extend the 
above Bloch sphere formalism to nanoelectrome- 
chanical systems. Our realisation is based on the 
two orthogonal fundamental flexural modes of a 
high quality factor nanostring resonator which are 
strongly coupled by a dielectric gradient field 9 . 
Full Bloch sphere control is demonstrated via 
Rabi 10 [ Ramsey 11 and Hahn echo 12 experiments. 
This allows manipulating the classical superpo- 
sition state of the coupled modes in amplitude 
and phase and enables deep insight into the deco- 
herence mechanisms of nanomechanical systems. 
We have determined the energy relaxation time 
Ti and phase relaxation times and T 2 % and 
find them all to be equal. This not only indi- 
cates that energy relaxation is the dominating 
source of decoherence, but also demonstrates that 
reversible dephasing processes are negligible in 
such collective mechanical modes. We thus con- 
clude that not only T\ but also T2 can be in- 
creased by engineering larger mechanical quality 
factors. After a series of ground-breaking exper- 
iments on ground state cooling and non-classical 
signatures of nanomechanical resonators in recent 
y ears [HHT3 this is of particular interest in the con- 
text of quantum information processing 1 18 em- 
ploying nanomechanical resonators 1 19 * 20 1 . 

While the dynamics of a two-level system under the 
influence of a pulsed external electromagnetic field was 
observed in atomic and nuclear spin physics decades ago, 
a mechanical analogon to such a system remained elu- 
sive for a long time. Only recently, coherent exchange 
of energy quanta between a mechanical and an electrical 
mode was achieved: In 2010, O'Connell et al.^ man- 
aged to control the swapping of a single quantum of en- 
ergy between a qubit and a mechanical resonator, while 
Palomaki et al.^ demonstrated the temporary storage 
of itinerant microwave photons in a mechanical resonator 



in 2012. At the same time, several approaches were em- 
ployed to achieve purely mechanical resonant coupling ei- 
ther between separat e reso nators^^ or different modes 
of the same resonator^H i n the classical regime. So far, 
the pulsed coherent control of the system was prevented 
by weak coupling, low quality factors or the lack of a 
sufficiently strong and fast tuning mechanism. 

We present the successful implementation of a purely 
mechanical two-level system with coherent time-domain 
control (see also the experiments independently per- 
formed at NTT using parametric coupling 25 ). To this 
end, we use a 250 nm wide and 100 nm thick, strongly 
stressed^ silicon nitride beam resonator with a length of 
50 um dielectrically coupled to a pair of electrodes used 
for detection^ as well as actuation and tuning^. The 
two fundamental flexural modes of the mechanical res- 
onator oscillating in the out-of-plane and in-plane direc- 
tion (see Fig. [TJ are coupled by cross-derivatives of the 
strong inhomogeneous electric field generated between 
the electrodes^. A constant dc voltage of -15 V is used 
to dielectrically tune the system close to the resulting 
avoided crossing, while the signals generated by an ar- 
bitrary waveform generator (AWG) enable time-resolved 
control vicinal to the anticrossing (see Fig. [T|d,c). Both 
voltages are added and combined with the rf actuation 
of the beam via a bias-tee and applied to one electrode. 
The other electrode is connected to a 3.6 GHz microstrip 
cavity, enabling heterodyne detection of the beam de- 
flection^ after addition of a microwave bypass capaci- 
tor at the first electrode^. These components as well 
as the mechanical resonator are placed in a vacuum of 
< 10 -4 mbar and cooled to 10.00 ± 0.02 K to improve 
the temperature stability as well as cavity quality factor. 
The microwave cavity is interfaced to the readout with a 
single coaxial cable and a circulator. 

When the system is driven by an external white noise 
source and the AWG output voltage is swept, the avoided 
crossing of the two modes shown in Fig.fTJ) can be mapped 
out, exhibiting a frequency splitting Q = 24, 249 ± 4 Hz. 
With a quality factor Q = -£j « 2- 10 5 and a linewidth of 
A/ « 40 Hz at the resonance frequency /, the system is 
clearly in the strong coupling regime of Af <C Q. For all 
measurements discussed in the following, an rf drive of 
-59dBm at 7.539 MHz, resonantly actuating the beam at 
an AWG voltage of V, is applied, which initialises the 
system in its in-plane mode (see black circle in Fig. [TJ^). 
A 1ms long, adiabatic voltage ramp up to 2.82 V brings 
the state to the point of minimal frequency splitting 
between the coupled modes. Here, the system dynamics 
is described by two hybrid modes formed by the in-phase 
and out-of-phase combinations of the fundamental flex- 
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FIG. 1. Nanoelectromechanical system, a, SEM mi- 
crograph showing oblique view of the 50 um long silicon nitride 
beam (green) and the adjacent, 1 um wide gold electrodes (yel- 
low), processed on top of the SiN. b, Electrical setup: the output 
of the arbitray waveform generator (AWG) and a dc tuning volt- 
age are added and combined with the rf drive via a bias-tee. The 
second capacitor acts as a bypass providing a uw ground path 
for the microwave detection^, which is connected to the other 
electrode, c, Resonance frequencies of the out-of-plane (out) 
and in-plane (in) mode of the resonator controlled by the AWG 
voltage at a constant dc tuning voltage of -15 V. The black circle 
marks the initialisation state at 0V and the frequency of the rf 
drive, while the green and blue circles correspond to the lower 
and upper state of the two-level system, respectively, separated 
by the frequency splitting 



ural modes. The adiabatic ramp thus transforms all the 
energy of the in-plane mode into the lower hybrid state, 
such that the two-level system, consisting of the two hy- 
brid modes, is prepared in its lower state. As the drive 
frequency remains constant (dashed line in Fig. the 
beam is no longer actuated and its energy is slowly de- 
caying. 

Now, the application of a continuous pum p tone with 
frequency Q will start Rabi oscillations^^ between the 
lower and upper state, as shown in Fig. [2] They can 
be measured directly by monitoring the time evolution 
of the output power spectrum at the frequency of one 
of the hybrid modes, here shown for the upper state at 
7.6028 MHz, and measured with a bandwidth of 10 kHz. 
All time-resolved measurements are averaged over 20 
(Rabi oscillations and Tl measurement) or 10 pulse se- 
quences (Ramsey fringes and Hahn echo). For a drive am- 
plitude of 100 mV (half peak-to-peak) we find a Rabi fre- 
quency of 8.3 kHz (see section II. A of the Supplementary 
Information for the frequency dependence of the Rabi os- 
cillations). In principle, the decay of these oscillations is 
governed by both energy relaxation, characterised by a 
rate 1/Ti, and phase decoherence, characterised by I/T2 
or 1/T 2 *, where T 2 * < T 2 includes reversible processes 
caused by slow fluctuations or spatial inhomogeneity of 
the coupling. For clarity, we use these well-known phe- 
nomenological constants in the same way as, e. g., in 
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FIG. 2. Rabi oscillations, a, Pulse scheme: the system 
is adiabatically tuned from the initialisation to the lower state, 
then a constant drive with frequency Q is turned on. b, The 
z projections of the decaying Rabi oscillations (data: dark blue; 
fit: red) can be directly measured with a spectrum analyser. The 
Bloch sphere in the inset shows the state of the Bloch vector at 
selected times, which are marked in the same colour in a. 
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FIG. 3. Energy relaxation, a, Pulse scheme: the system 
is adiabatically tuned from the initialisation to the lower state. 
An additional 7r-pulse is used to rotate it to the upper state, b, 
Measured exponential decay of the lower (data: green; fit: dark 
green) and upper (data: blue; fit: dark blue) state. The Bloch 
sphere in the inset shows the state of the Bloch vector at selected 
times, which are marked in the same colour in a. 



spin systems^, as discussed in more detail in the Supple- 
mentary Information section I. 

The exponential decay of a state's energy defines T\. 
The corresponding measurement is shown in Fig. [3] for 
both lower and upper state: The system is once again 
prepared in the lower hybrid state. To reach the upper 
state, a subsequent 7r-pulse is applied, thus performing 
one half of a Rabi cycle which transfers the system to 
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the upper state (see Supplementary Information section 
II for details on the frequency and amplitude calibra- 
tion of the applied pulses). The exponential decay is 
then measured directly with a spectrum analyser at a 
bandwidth of 3 kHz, exhibiting different relaxation times 
T 1: i = 4.83 ±0.1 ms and T ljU = 4.02 ±0.1 ms for the lower 
and upper mode, respectively. They correspond to the 
spectrally measured quality factors. Previously, it has 
been shown that, at maximum coupling, the two hybrid 
modes should have the same quality factor and thus T\ 
timeBD. However, both modes are affected by dielectric 
damping 28 , leading to the observed difference. 

To measure the T 2 * time, a 7r/2-pulse is used after the 
preparation in the lower state to bring the system into 
a superposition state between lower and upper hybrid 
mode. The frequency of the pulse is detuned to Q ± 
500 Hz, leading to a slow precession of th e state vector 
around the z-axis of the Bloch sphere^^. As a result, a 
second 7r/2-pulse after time r does not always bring the 
system into the upper state, but a slow oscillation, the 
so-called Ramsey fringes, is observed when the delay r 
between the two pulses is varied and the z-projection of 
the state vector is measured after the second pulse, as 
shown in Fig. [4j The decay constant of this oscillation is 
T 2 *, while the decay of the mean value corresponds to an 
effective T\ of both modes. The fit in Fig. [4]d results in 
T 2 * = 4.44 ± 0.1 ms and T x = 4.31 ±0.1 ms. The energy 
relaxation time of the superposition state T x is identical 
to the reciprocal rate average of the two hybrid modes 
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as the mechanical energy oscillates between the two 
modes with frequency Q (see Supplemental Video). 

By including an additional 7r-pulse at r/2 into the 
Ramsey pulse scheme and replacing the final 7r/2-pulse 
by an 37r/2-pulse to once again rotate to the upper state 
(see Fig.pl), the T2 time can be measured in a Hahn echo 
experimental^. The 180° rotation flips the state vector 
in the xy-plane of the Bloch sphere, thus reversing the ef- 
fects of a fluctuating or inhomogeneous coupling strength 
Q in the second delay interval of r/2 and thereby cancel- 
ing their contribution. The frequency of the pulses is once 
again exactly £7, as all three pulses need to be applied ex- 
actly around the same axis. The resulting decay curve 
represents T2, for which a value of T2 = 4.35 ±0.1 ms can 
be extracted from the fit in Fig. [5Jd. 

The good agreement between T2 and T 2 * clearly shows 
that reversible elastic dephasing, e. g. caused by tem- 
poral and spatial enviromental fluctuations or spatial in- 
homogeneities, does not noticeably increase decoherence. 
Although the experiment is performed with billions of 
phonons, they all reside in the same collective mechanical 
mode and thus all experience an identical environment. 
This strongly constrasts the behaviour found e. g. in 
spin qubits, where the hyperfine interaction with ener- 
getically degenerate nuclear spins causes T 2 * <C T2 (see 
Supplemental Information section III). 
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FIG. 4. Ramsey fringes, a, Pulse scheme: the system is 
adiabatically tuned from the initialisation to the lower state. A 
7r/2-pulse creates a superposition state, and after a delay r a sec- 
ond 7r/2-pulse is applied, b, A 500 Hz detuning between drive 
and precession frequency leads to a slow rotation of the super- 
postion state in the equator plane of the Bloch sphere, giving 
rise to a beating pattern in the measured z component after the 
second pulse (data: dark blue; fit: red). The Bloch sphere in 
the inset shows the state of the Bloch vector at selected times, 
which are marked in the same colour in a. 
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FIG. 5. Hahn echo, a, Pulse scheme: the system is adi- 
abatically tuned from the initialisation to the lower state. A 
7r/2-pulse creates a superposition state, and after a delay r/2 a 
7r-pulse mirrors the state vector to the other half of the Bloch 
sphere. After another delay of r/2 a 37r/2-pulse is used to rotate 
to the upper state, b, The inverse evolution of the system during 
the two delay times cancels out any broadening or slow preces- 
sion effects, thus the system always ends up along the z-axis 
and no oscillation is observed (data: dark blue; fit: red). The 
Bloch sphere in the inset shows the state of the Bloch vector at 
selected times, which are marked in the same colour in a. 
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It is more surprising that the phase coherence time T2 
is equal to the average energy relaxation time T\. This 
indicates the absence of measurable elastic phase relax- 
ation processes in the nanomechanical system, such that 
the observed loss of coherence is essentially caused by 
the energy decay of the mechanical oscillation (see also 
Supplemental Information section III). Earlier research^ 
suggests that the dominant relaxation mechanism in sil- 
icon nitride strings is mediated by localised defect states 
of the amorphous resonator material, described as two- 
level systems at low temperature. They facilitate en- 
ergy relaxation by providing the momentum required to 
transform a resonator phonon into a bulk phonon. For 
this process to lead to elastic phase relaxation, an ex- 
cited defect state would have to re-emit the phonon back 
into the resonator mode, which is extremely unlikely due 
to the weak coupling between the two. In conclusion, we 
demonstrate coherent electrical control of a strongly cou- 
pled (CI ^> ^) nanomechanical two- level system, employ- 
ing the pulse techniques well-known from coherent spin 
dynamics in the field of nanomechanics. Each superposi- 
tion state of the two hybrid modes on the Bloch sphere 
can be addressed by a sequence of the described pulses. 
The presented system stands out by the finding that the 
elastic phase relaxation rate 1^ is negligible compared to 
the energy decay rate , leaving room for improvement 
of the coherence via increased quality factors. 

In light of the recent breakthrough in ground-state 
cooling of nanomechanical resonators'^"^, the coherent 
manipulation schemes presented here open new applica- 
tions for nanomechanical systems in quantum informa- 
tion. Not only can they be used as efficient interfac es for 
quantum state transfers in hybrid quantum systems^^^, 
but by creating coupled, quantised resonators^ quantum 
computations can be carried out directly using nanoelec- 
tromechanical two-level systems^. 
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SUPPLEMENTAL MATERIAL TO "COHERENT 
CONTROL OF A NANOMECHANICAL 
TWO-LEVEL SYSTEM" 

I. RELAXATION TIMES 

In the main text, the relaxation constants T\ and T2 
are introduced phenomenologically to define the expo- 
nential decay times extracted from the energy relaxation 
and Hahn echo experiments. But different from the two- 
state spin systems associated with experiments on the 
Bloch sphere, the mechanical system investigated here 
actually has three states: the two coupled hybrid modes, 
i. e. the lower and upper state as well as the ground 
state (the thermally occupied phonon bath) into which a 
phonon can relax from either mode, see Fig. [Sl| 



upper 



lower 



bath 



FIG. SI. Levels of the mechanical system: Schematic 
representation of the phonon bath and the upper and lower 
hybrid state of the mechanical system, separated by the fre- 
quency splitting Q, with respective decay constants k u and 

For one, this makes it necessary to prepare the lower 
state prior to any measurement, as none of the states 
of the nanomechanical two-level system is automatically 
populated (except for the comparatively weak thermal 
excitation). Furthermore, it introduces additional terms 
to the Bloch equations: Assuming two independent decay 
rates k u and k\ for the upper and lower state and follow- 
ing reference 31 , the rotating-frame Bloch equations can 
be written as 



M x (t) = 



My(t) 



M z (t) 



1 



k)M x (t) + AM y (t) 



(SI) 



i/t 2 



Tj_ 



k M y (t) - AM x (t) - uj R M z (t) 



i/t 2 



k)M Z (t)+LJ R My(t) 



(S2) 
(S3) 



l/Ti 



Here, A is the detuning between drive frequency and 
coupling strength Q and uor reflects the drive strength 
and corresponds to the frequency of Rabi oscillations. 
k = ku + kl is the average decay rate, T± the relaxation 
time in the equator plane of the Bloch sphere, Tji the 
relaxation time along the z direction and M x , M y and 




0.0 

in-plane displacement 

FIG. S2. Oscillation in the superposition state: The 

normalized displacement of the resonator is calculated for one 
full precession period in the superposition state, i. e. for 
equal amplitudes of the lower and upper state. This phase 
trajectory is plotted in the basis of the in-plane and out-of- 
plane mode (horizontal and vertical axis of the plot). 



M z denote the respective components of the state vec- 
tor. The phenomenological parameters l/Ti and I/T2 
can be identified as the sum of the respective rates. 

The measured values of Ti are consistent with the qual- 
ity factors of the corresponding modes, which are only 
limited by k u and fcj. This implies a negligibly small 
1/Ty. As the measurements show that T\ is equal to 
T2, 1/Tj_ must also be negligible. Thus, the two coher- 
ence times Ti and T2 are solely limited by the average 
mechanical damping of the two resonator modes. 

Reference 31 also introduces an additional term in 



equation (S3) taking into account the difference between 



the two decay rates k u and k\. It leads to a slow tilt of 
the state vector of a superposition state away from the 
equator plane towards the state with the smaller decay 
rate. However, this effect plays no role in the classical 
system presented here: In a classical superposition state, 
the energy is distributed between both modes. The sys- 
tem performs an oscillation between its two fundamental 
modes, changing the direction of rotation with the pre- 
cession frequency. Thus, the time spent in each mode is 
equal, and the state just experiences the average decay 
constant. A plot of the phase space trajectory for one 
precession cycle is shown in Fig. |S2[ using a 10 times ex- 
aggerated coupling strength. The horizontal and vertical 
axes correspond to an in-plane and out-of-plane oscilla- 
tion, while a diagonal motion is associated with the lower 
and upper hybrid mode. An animated version of this plot 
is available as an ancillary file on arXiv.org. 
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II. PULSE CALIBRATION 



B. Pulse length and spacing 



The first step in characterising the system is to mea- 
sure an avoided crossing as shown in Figure 1 of the main 
text. From its fit, the approximate frequency splitting Q 
and the AWG voltage required to adiabatically tune to 
the lower state is extracted 9 . 



A. Pulse frequency 

To precisely determine the correct pump tone, the fit- 
ted frequency splitting is not accurate enough. Instead, 
the frequency of Rabi oscillations is monitored while 
sweeping the pump frequency. The quadratic dependence 
for small detunings 2 , as shown in Fig. S3 1 allows to fit the 
measured points and extract the lowest Rabi frequency 
and thus the exact pump frequency corresponding to zero 
detuning. 
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FIG. S3. Pump frequency tuning: The frequency of Rabi 
oscillations depends quadratically on the pump frequency de- 
tuning. The upper plot displays the measured signal power 
colour- coded versus pump frequency and measurement time. 
The extracted Rabi frequencies (blue points: data; red line: 
parabolic fit) are shown in the lower plot. The minimum of the 
parabola corresponds to the zero-detuning pump frequency. 



The pulses applied to the system should be as short 
as possible to allow fast control sequences. The lower 
limit of the pulse length is one period of the pump signal, 
as abrupt voltage jumps of a chopped up sine wave will 
disturb the system. The shortest applied pulse (a tt/2 
pulse) is thus set to a duration of 1/ (pump frequency) 
so that it consists of exactly one sine wave. The other 
pulses are correspondingly longer, as shown in the pulse 
schemes in Figures 3-5 of the main text. 

The delay time between two pulses in Figure 4 and 
5 also needs to be a multiple of 1/ (pump frequency) , 
as otherwise the second rotation will not be carried out 
along the same axis. This is apparent from the Ramsey 
fringe experiment, where the observed beating pattern, 
and thus the effective rotation angle, is caused by an 
intentional detuning of 500 Hz. 



C. Pulse amplitude 

After both length and frequency of a pulse are fixed 
as described above, its amplitude has to be adjusted to 
define the rotation angle achieved with each pulse. To 
this end, four sine periods (i. e. a 4 • tt/2 = 2tt pulse at 
the desired amplitude) are applied to the system, vary- 
ing the amplitude. The population of the upper state 
is measured. At zero amplitude, the pulse has no effect 
and all energy remains in the lower state. With increas- 
ing amplitude, the achieved rotation angle increases, and 
the first minimum corresponds to the desired 2tt charac- 
ter of the pulse. The population P u of the upper state 
can be described as 



P u = sin — — 
An 



(S4) 



where A is the pulse amplitude and Aq the amplitude 
corresponding to a rotation of 2ir. A fit to the measured 




100 



200 300 
Pulse amplitude (mV) 



400 



FIG. S4. Pulse amplitude sweep: The population of the 
upper state after quadruple- sine- wave pulses of different am- 
plitudes (blue points: data; red line: fit) demonstrates the 
correct behaviour for an amplitude of 231 mV. 
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data, as shown in Fig. [S4j can be used to extract the 
correct pulse amplitude. This measurement has to be 
repeated with the detuned pump frequency used for the 
Ramsey fringes experiment, as the slightly different fre- 
quency also leads to a small shift in the amplitude of the 
pulse. 



III. PHASE RELAXATION 

To define the elastic phase relaxation rate two situ- 
ations have to be distinguished: In systems where energy 
relaxation occurs only from the upper to the lower level 
(e.g. references 8 ' 32 ) the phase relaxation rate is defined 



D 31,33,34 



1 



1 



T 2 2Ti " 



(S5) 



In contrast to that, in a system dominated by sponta- 
neous decay to a third state and no relaxation between 
the upper and lower level 33 ' 34 (i. e. 1/Tn = and fi- 
nite fc, see equation (S3)), the phase relaxation rate T^, 
in this case equal to ljT± (defined in equation (SI) and 
equation (|S2])), is given as 



1 1 



(S6) 



This is the case for the system presented here. As T\ 
and T2 are equal within the measurement acurracy, 
can not be determined from the experiment. This shows 
that the measured phase decoherence is solely caused by 
energy relaxation. Processes changing the phase but pre- 
serving the state's energy seem to play no role. 

To compare this mechanical two-level system to other 
coherent systems, it helps to take a general look at the 
possible decay processes: Inelastic processes in which en- 
ergy is transfered to a thermal bath are irreversible. They 
are directly represented in the T\ time and also pose a 
limit to T2 via the two above equations S5 or S6 Ir- 



reversible elastic interactions lead to a non-zero and 
thus reduce T2, whereas reversible phase decay processes 
can be measured and controlled e. g. by a Hahn echo 
experiment and only decrease T 2 *. 

In most coherent nanoscale solid-state systems, the 
coupling to a fluctuating thermal bath of phonons, pho- 
tons, two- level systems or (nuclear) spins leads to one 
or more of the above processes. For example, supercon- 
ducting qubits 35-37 suffer from flux, charge and pho- 
ton noise. In gate-defined spin qubits e. g. nuclear 
spin 38-40 and phonon 41 interactions limit the perfor- 
mance, whereas NV centres in diamond 42 couple to the 
surrounding nuclear spin bath. In the amorphous dielec- 
tric system presented here, relevant loss mechanisms oc- 
cur via defect states with a broad energy spectrum often 
associated with two-level systems 43 , and the phononic 
environment, as long as no additional electronic noise is 
introduced via the measurement devices and tuning volt- 
ages. As the mechanical modes under investigation are 



situated within the suspended beam, they effectively re- 
side inside a phonon cavity and couple extremely weakly 
to the phonon bath of the bulk sample via the narrow 
clamping points 44 . The exchange of energy of the dis- 
crete long-wavelength resonator modes and the conti- 
nous shorter-wavelength phonon spectrum of the beam 
is found to be mediated by the defects. As only higher 
energy phonons with small wavelengths can effectively 
transmit energy through the clamps to the bath, scatter- 
ing of thermally excited higher energy phonons at defect 
states in a three particle interaction is the most likely 
process. These inelastic processes destroy energy as well 
as phase and likely explain why we find T\ = T^. 
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